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Abstract 



The distributions of the angular transmission coefficient and of the total 
transmission are calculated for multiple scattered waves. The calculation is 
based on a mapping to the distribution of eigenvalues of the transmission 
matrix. The distributions depend on the profile of the incoming beam. The 
distribution function of the angular transmission has a stretched exponen- 
tial decay. The total-transmission distribution grows log-normally whereas it 
^ | decays exponentially. 

42.25-p,78.20.Dj,72.10.-d,05.40+j 

Transport through mesoscopic systems is a wide field of interest and is studied with 
a variety of waves: sound, microwaves, electrons, and light. The wave character leads to 
interference between the transmission channels. This causes large fluctuations. Therefore 
the observables are not always characterized by mean values, but their entire distribution 
functions are of interest. This is particularly prominent in the distribution of eigenvalues of 
the transmission matrix. The eigenvalues have a bimodal distribution peaked around and 
1 with a mean value £/L <C 1, where I denotes the mean free path and L the thickness of 

medium the transmitted signal consists of speckles. These speckles are the first quantity 
of interest. The intensities of the speckles obey in first approximation the Rayleigh law, 
which is a negative exponential distribution. In the mesoscopic regime interference modifies 
this distribution. The leading correction was derived by Shnerb and Kaveh [|J] . Genack and 
Garcia observed a deviation from the Rayleigh law at large intensities ||. A crossover to 
stretched exponential behavior was derived by Kogan et al. ||. 

The second quantity is the total transmission. It is obtained by integrating over the 
outgoing surface, thus collecting all speckles. When adding the speckles of the N independent 
channels, the law of large numbers predicts a Gaussian distribution with variance of the order 
1/N. However, interference broadens the variance to L/(£N) 0. This was confirmed in a 
recent experiment of de Boer et al. ||. Also in computer simulations by Edrei et al. the 
Gaussian behavior was observed, while a cross over to a log-normal distribution was seen for 
large disorder [[J . De Boer et al. also measured a small but clear deviation from a Gaussian, 
due interference of three transmission channels. 



the sample 

Three distinct transmission quantities can be measured. If a laser illuminates a diffuse 
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The third quantity which can be extracted is the conductance, which is obtained by sum- 
ming also over all incoming directions. An extensive discussion of the distribution function 



of the conductance including a review was given by Altshuler et al. [TO 



These three quantities are dominated by different interference processes. This can al- 



ready be seen by considering their autocorrelation functions ]II,I^fl . Let T a t be the angular 
transmission coefficient for a speckle spot b arising from incoming wave a. Its correlator is 
called the C\ correlation function. It is dominated by a disconnected diagram. The corre- 
lator of the total transmission T a = J2bT a b is called the C2 correlation, and dominated by 
the interference of two diffusons. In a diagrammatic approach it is a connected, loopless dia- 
gram. The interference process mixes two incoming into two outgoing diffusons. Finally, the 
correlator of the conductivity T = Ylab T a b is called C3 correlation or Universal Conductance 
Fluctuation. It is dominated by a connected loop diagram. 

We shall first calculate the distribution function of the total transmission. Subsequently 
we apply these ideas to the statistics of the angular transmission coefficient. 

Consider a three dimensional slab with dimensions W x W x L, (W ^> L) with elastic 
scatterers at quenched random positions. We consider the 'metallic' regime, where the 
dimensionless conductance g is large. Let us look at the transmission problem for a scalar 
plane wave with unit incoming flux in channel a, 

CX r ) = ~J%£= eX P( ic laP + ikfl a z) Z < 0, (1) 

where A = W x W denotes the area of the slab, k is the wavenumber, p = (x, y) is the 



transversal coordinate and ji a = Jl — q 2 /k 2 = cos# a , where 6 a is the angle with respect 



to the z— axis. The precise form of the propagators was derived in Ref. ||13|| . An incoming 
diffuson has the form 

r in,s _ 47T7i(/A a ) L-Z 

LAT) ~ ktAp a L ' {Z) 

where T\ describes the limit intensity of a semi-infinite system. We shall term the propagator 
summed over all incoming channels the 'total-flux diffuson': 

nr),ECW = y^ ( 3 ) 

In the bulk the difference between a total-flux diffuson and a diffuson arising from plane 
wave incidence is just described by an overall factor e a , 

C(r)=6 a £ in (r), (4) 

with e a = 7"i(/i a )7r / '{k 2 Afi a ) satisfying the sum rule Y,a e a = 1- Similarly, the outgoing 
diffusons are 

Cr(r) = e,£ out (r), = £ £. (5) 

The wave is transmitted into outgoing channel b with transmission amplitude t a b = 
2k^/p a fibG a b an d transmission probability T ab = \t a b\ 2 . The average total transmission is 
obtained by summing all outgoing channels 
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(T a ) = (T,T ab ) = T -^. (6) 



The average conductance is given by 



h 2 Af 



thus (T a ) = e a g, while one also has (T ab ) = e a e b g [JT3[ . 

We consider the j th cumulant of T a . In a diagrammatic approach this object has j 
transmission amplitudes t ab and an equal number of Hermitean conjugates t\ a = t* ab . Ex- 
plicit calculation for the second and third cumulant showed that the leading diagrams are 
connected, but have no loops ||. We assume that this holds to any order. Let us fix the 
external diffusons in the term t abl t\ ia t ab2 ■ ■ •i<%4 a - Contributions to the sum over hi only 
come from diagrams with outgoing diffusons that have no transversal momentum. These are 
the diagrams where the lines with equal 6j are paired into diffusons. The outgoing diffusons 
are fixed now. For connected loopless diagrams there are (j — 1)! distinct choices for pairing 
the incoming amplitudes. Next we factor out the incoming and outgoing diffusons and group 
the remainder of the diagrams into a skeleton K. Making use of Eq. @ we obtain: 

ra c on=4(j-l)! / dndri ■ • •dr,dr;./: in (r 1 )/: out (r' 1 ) 

.•.^(r^^O^r^r;,-..,^^.). (8) 

The integral just describes 

(Tr(ttt)i) = £ (t aib A im t a2b2 ■ --t^tlJ. (9) 

ai, bi,—, a j,bj 

Indeed, for this quantity there is only one way to attach incoming and outgoing diffusons to 
K. The sums over the indices lead exactly to the total- flux diffusons in Eq. We thus 
find 

TOcon = (j - l)!e{(Tr(ttt)i), (10) 

which is the crucial step in the derivation. The eigenvalues T n of the transmission matrix vt 
can be expressed as T n = 1/ cosh 2 L7„. Under very general conditions [jl4[ the distribution 
of the Lyapunov coefficients 7 n is uniform This implies that: 

Normalizing with respect to the average, we introduce s a = T a / (T a ). The generating function 
of the connected diagrams is easily calculated 

$con(2) = E ' ~ ( S a)con 

3=1 J' 



g\og 2 Ul+x/g + ^r g ). (12) 
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Since the cumulants are solely given by connected diagrams, the distribution of s a follows 

as 

/ioo oj<£ 
— exp [xs a - $ con (^)] • (13) 
-ioo Zlll 

For s a near unity and large g we can expand $ up to order x 2 , recovering the Gaussian 
behavior found by Kogan et al. ||: 

p(sa) * \ZS exph f {sa - i)2] (i4) 

The integrand in Eq. (|H|) has a branch cut from x = —g to x = — oo. For s a < the 
contour can be closed to the right and p(s a ) vanishes. The shape for small s a (and large g) 
is dominated by a saddle point. One finds essentially a log- normal growth: 



p(s a ) ~ exp 



g g / 2 , , 2 
4~ 4 l 1 ° S ^ + l0Sl0g ^- 1 



(15) 



Also for large s a we can apply steepest descent. Here one finds a simple exponential decay 

7T 2 

p(s a ) ^ exp(-gs a + g—), s a > 1. (16) 

In Fig. m we present the distribution ( |13D for some values of g. At moderate g the deviation 
from a Gaussian is clearly seen. 

So far we have considered the case of an incoming plane wave. In optical systems a Gaus- 
sian intensity profile is more realistic. For perpendicular incidence the incoming amplitude 
is -^ in (r) = W~ l J2 a 0(9a)V'a n ( r )) where VC( r ) i s the plane wave of Eq. ([I]), and where 

<f>(q a ) = V2^p exp(~p 2 Q q 2 a ). (17) 

We consider the limit where the beam is much broader than the sample thickness (po ^ L) 
but still much smaller than the transversal size of the slab (po W). (A smaller beam 
diameter complicates the problem; when incoming transverse momenta, which are of order 
1/po, become of the order of 1/L, the diffusons will take a different form |5J). Here the 
momentum dependence of the diffusons can be neglected. Due to integration over the center 
of gravity, each diagram involves a factor Ad^q^q 1 - In the j th order term there occurs a 
factor 

F i = % E 0(gi)0*(gi)---0fe)0*(^)^ 



d 2 p|0(p)| 2j - (18) 
For a plane wave we have |0(p)| = yA, and Fj = A X ~K For our Gaussian beam we obtain 
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It is thus convenient to identify Aq = \np\ with the effective area of a Gaussian beam. As 
compared to the plane wave case, the j th order term is smaller by a factor 1/j for a Gaussian 
profile. This implies for the generating function of the connected diagrams 



•-<*> = 'j!,7 v lf + 7 + V,J- (20) 

For small s a (and large g) there is again a log-normal saddle point. For large s a the dominant 
shape of the decay is given by the singularity at x — —g and again yields p(s a ) ~ exp(— gs a ). 
In Fig. we present the distribution function for different values of g. 

By expansion of the generating function we recover the results for the second and third 
cumulant, obtained previously by de Boer et al. || 

(■S a )cum — g~) (s o )cum = y( S n)cum (21) 

for a plane wave. For a Gaussian beam we recover 

/ 2\ / 3\ 16 J 2\2 /oo\ 

\ s a /cum = ; \ s a /cum = 7~\ S a/cuni - \^) 

3g 5 

These last relations were confirmed experimentally M. 

We apply the same method for the distribution of the angular transmission coefficient. 
In the plane wave situation the average reads (T^) = e^g. Let us count the number 
of connected loopless diagrams that contribute to T ] ah = t ao t\ a t a b ■ ■ -t ab tl a . In this case all 
pairings into outgoing diffusons contribute. This yields an extra combinatorical factor j\ in 
the j th moment: 

(Ti) CO n=jKj-iy44m j t) j ). (23) 

For the normalized angular transmission coefficient s a t, = T a b/(T a fc), we introduce the fol- 
lowing generating function of the connected diagrams: 

r-r- 



*con(a;) = $con(a;), (25) 



It is easy to see that 



with $ con given by Eq. fll^D for plane wave incidence and by Eq. (pD| ) for a broad Gaussian 
beam, respectively. In contrast to the total-transmission distribution, the cumulants are not 
only given by the connected diagrams. Kogan et al. showed that the summation of the 
disconnected diagrams can be done elegantly by performing an additional integral. Using 
this result one gets 

. . f 00 dv f io ° dx ( , A , , 

P{Sab) = — / — exp h xv - W con (x) . (26) 

Jo v J~ioo 2m V v J 
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The speckle intensity distribution is plotted in Fig. for an incoming plane wave. 

For large g and moderate s ab we have $ con (x) ~ x and we recover the Rayleigh law: 
p(sab) = exp(— s a b). The leading correction is found by expanding in 1/g 



]>{■%>>) - " 



1 + ^(4,-4^ + 2) 



(27) 



This was derived previously by Shnerb and Kaveh ; here we have related the prefactor of 
the correction term to the conductivity. Genack and Garcia fitted their data to this relation 
and found g = 14.6. Our Eq. ( p6|) describes these data very well for g = 14.4. 
For large s a b one can again apply steepest descent, which yields 

p(s ab ) ~ exp(-2y/gs ab ). (28) 

This stretched exponential tail differs from the form p(s ao ) ~ exp[— (&lgs 2 ab /lQ) l l' i ] asserted 
by Kogan et al. Their findings are based on truncating $ con (x) after order x 2 , which corre- 
sponds to including only the simplest connected diagram. Taking the full generating function 
into account, we find a qualitatively different saddle point. 

A Gaussian profile of the incoming beam leads to a different distribution with the same 
asymptotic behavior. 



Our results remain valid if there is a mismatch in the indices of refraction |R| , but may 
change if the system is optically not very thick. 

In conclusion, we have calculated the distribution function of the total amount of light 
transmitted through a disordered slab. It was related to the distribution function of the 
eigenvalues of the transmission matrix, which is known since the Lyapunov coefficients are 
uniformly distributed. Comparing to previous diagrammatic calculations M it became clear 
that this uniform distribution comes from summing connected, loopless diagrams. 

Next we have calculated the distribution of the angular transmission intensity. Here a 
deviation from Rayleigh statistics is found, of new stretched exponential form. Our calcula- 
tions are in agreement with measurements of both total and angular resolved transmission. 
As our results are based on loopless diagrams, they constitute the mean field expressions. 
Near the Anderson transition loop diagrams will change the distribution functions. 
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FIGURES 



FIG. 1. Intensity distribution of the total transmission, in units of yj3g/4n, versus the normal- 
ized intensity s a for an incoming plane wave. With g = 2,4, 16, and 64 (upper to lower curve). 

FIG. 2. Intensity distribution of the total transmission, in units of ^figj 7 ^, versus the normal- 
ized intensity s a for an incoming wave with Gaussian profile. With g = 2,4, 16, and 64 (upper to 
lower curve). 

FIG. 3. Intensity distribution of speckles versus the normalized intensity s a b for an incoming 
plane wave. With g = 2,4, and 8 (upper to lower curve at s a b = 5). The dashed line corresponds 
to the Rayleigh law (g = oo). 
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